Correlation of inflation-produced magnetic flelds with scalar fluctuations 



Robert R. Caldwell and Leonardo MotteQ 
Department of Physics & Astronomy, HB 6127 Wilder Lab, Dartmouth College, Hanover, NH 03755, USA 

Marc Kamionkowski 

Calif orrna Institute of Technology, Mail Code 350-17, Pasadena, California 91125, USA and 
Dept. of Physics & Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA 

(Dated: September 22, 2011) 

If the conformal invariance of electromagnetism is broken during inflation, then primordial mag- 
netic fields may be produced. If this symmetry breaking is generated by the coupling between 
electromagnetism and a scalar field — e.g. the infiaton, curvaton, or the Ricci scalar — then these 
magnetic fields may be correlated with primordial density perturbations, opening a new window 
to the study of non-gaussianity in cosmology. In order to illustrate, we couple electromagnetism 
to an auxiliary scalar field in a de Sitter background. We calculate the power spectra for scalar- 
field perturbations and magnetic fields, showing how a scale-free magnetic field spectrum with rms 
amplitude of ~ nG at Mpc scales may be achieved. We explore the Fourier-space dependence of 
the cross-correlation between the scalar field and magnetic fields, showing that the dimensionless 
amplitude, measured in units of the power spectra, can grow as large as ~ bOOHj /M , where Hj is 
the inflationary Hubble constant and AI is the effective mass scale of the coupling. 



I. INTRODUCTION 

The predictions of the simplest single-field slow-roll 
models of inflation agree remarkably well with current 
cosmological data, yet experience gained from effective 
field theories suggests that this model is likely not the 
whole story. A vast literature has now arisen to explore 
ultraviolet completions and their predictions for future, 
more sensitive, observations [HE]- One of the principle 
lines of investigation has been the predictions for non- 
gaussianity due to self-couplings, nontrivial infiaton ki- 
netic terms or interactions between multiple fields asso- 
ciated with inflation [SHS]. 

Another possibility for beyond single field slow roll 
physics is coupling of the infiaton, or some other specta- 
tor field, to electromagnetism. If such a coupling breaks 
the conformal invariance of electromagnetism, then quan- 
tum fluctuations in the electromagnetic field may be am- 
plified into classical magnetic fields in much the same 
way as quantum flucutations in the infiaton (graviton) 
become density perturbations (gravitational waves). It 
has been suggested that such inflation-produced mag- 
netic fields may provide the seed fields required for galac- 
tic dynamos [6l-fT4] , but it may also be that the signatures 
of such magnetic fields may be observed in the cosmic mi- 
crowave background [TSHSH], and thus shed light on in- 
flation, even if they are unrelated to galactic magnetism. 
Either way, the search for primordial magnetic fields pro- 
vides an additional observational probe of the physics of 
inflation to parallel that obtained from non-gaussianity 
searches. 

Here we explore the cross-correlation between primor- 
dial magnetic fields and a scalar field in a toy model 
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in which the scalar field is coupled to electromagnetism, 
with no gravity, in a fixed de Sitter background. The ho- 
mogeneous time evolution of the scalar field breaks the 
conformal invariance of electromagnetism. We first cal- 
culate the quantum mechanical spectrum of scalar- and 
magnetic-field fluctuations produced, and we then cal- 
culate the cross-correlation between the scalar and mag- 
netic fields. 

If the scalar field is a curvaton field, and if that cur- 
vaton is responsible for primordial perturbations, then 
the scalar-field-magnetic-field cross correlation we cal- 
culate will be precisely the density-magnetic-field corre- 
lation observed in the Universe today. Our calculation 
also illustrates the principal ingredients that will arise in 
a density-pcrturbation-magnetic-field correlation if the 
scalar field is the infiaton. 

In Section |IT] we introduce our model, work out the 
dynamical behavior, and evaluate the two-point statis- 
tics of the scalar and magnetic fields. In Section III 



we 



present the calculation of the cross-correlation, and we 
analyze its behavior in Sec. |IV| We conclude in Sec. |V] 
Throughout, we work in spatially-flat Robertson- Walker 
coordinates, with line-element ds^ = a'^{ri){—dr]'^ + djp). 



II. MECHANISM OF MAGNETIC FIELD 
AMPLIFICATION 



The action for our model is 
S= I d*x 



where (j){x, t) is the scalar field, and F^j^ the electro- 
magnetic field-strength tensor. The scalar-field poten- 
tial is V{4>) = —3nMHj(j), and the coupling function is 



1^(0) 



We suppose that some other field is 



2 



driving inflation. In practice we consider a fixed de Sit- 
ter background with Hubble constant Hi, whereby the 
scale factor is a{rf) = ~l/{riHj) for the run of conformal 
times —oo<r]<r]i<0 and 77/ marks the end of infla- 
tion. Aspects of this model have previously been studied 
[21 [12] , but we revisit the details in preparation for our 
later calculations. 



Scalar Field 



The scalar-field equation of motion is 



□0 = 



dV 1 dW 



and it has a solution, 

= Co + ciif - nMln{r]/rii), 



(1) 



(2) 



where we assume there is no homogeneous electric or 
magnetic field. We take the integration constants cq 
and ci to vanish so that W{(j)) = 1 at the end of infla- 
tion. In this way, the usual electromagnetic Lagrangian 
is recovered for the post-inflationary epoch, and we as- 
sume that some mechanism stops the subsequent evolu- 
tion of (j), so that the standard Maxwell equations are 
preserved at all times after inflation. We also define 
I{r]) = [W {(j}{r]))]^^^ — [rj/rji)^"', which will appear in 
our analysis of the electromagnetic field. 

The scalar field has fluctuations (5(/)(x, 77), about its ho- 
mogeneous component, described by the evolution equa- 
tion. 



+ 2nS(f,' + {a^V^^ - \7^)S(f, = 0, 



(3) 



where H — a' /a and is the spatial Laplace operator. 
To be clear, we fix the background to be pure de Sitter 
spacetime and subsequently "turn off" gravity, so that 
there are no fluctuations of the spacetime metric. Fol- 
lowing standard procedures, the quantized scalar field is 
decomposed in terms of time-dependent mode functions 



d(l){x,r]) 



(2^ 



(4) 



where and a}^ are respectively annihilation and cre- 
ation operators that satisfy [afe,a|,,] — {2TT)^S{k — k'). 
The uncertainty relation for the scalar field and its con- 
jugate momentum 54>' , 



[S(j){x,T]),S(j>'{y,ri)] = i5{x-y)/a^{'q), 



(5) 



results in a constraint to the two linearly independent 
solutions to the mode equation. Because the effective 
mass is zero, Vrf,^ = 0, we obtain the solution. 



Hi i — krj 



-ikr) 



(6) 



corresponding to the Bunch-Davies state, having positive 
frequency in the remote past, 77 — )■ —00 for fc|?7| ^ 1. 
The requirement 5(j> <^ 4> that the fluctuations are small 
translates into the bound Hi/M <^ 1. Finally, the two- 
point correlation function is 



{S(t){x,ri)S(l){y,r])) 



where the scalar-field 



d^k 
(27r)3 



p^^(^) (7) 



power spectrum — defined by 
- k')Ps^{k) and <5d is the Dirac 
is Ps(p{k) = Hj/2k^, valid for modes out- 



{27rf6D{k 



delta function- 
side the horizon at the end of inflation. 

The root-mean-squared amplitude — the correlation 
function at zero lag (at x = y) — is divergent at both 
the infrared and ultraviolet limits. Hence, we bound the 
run of wavenumbers to [fcmin, fcmax], so that 

= ((<50)')'^' = ^ (lnA;,„ax/fcmin)^ (8) 



gives the rms scalar-field fluctuation. In practice, we as- 
sociate the minimum wavenumber with the present-day 
Hubble radius — i.e., fcmin = 2t:Hq — and the maximum 
wavenumber with an astrophysical scale that we indicate 
by A. 



B. Electromagnetism 

The full action for electromagnetism includes not only 
the free Maxwell field, but also the coupling to charged 
particles as well as the action for the charged particles 
themselves. Including these additional terms, we may 
write 



X V-5 



where Cq is the Lagrangian for charged particles. 



(9) 
The 



electromagnetic coupling, or electric charge, is inversely 
proportional to I{ri). Consequently, in the case n > 
the coupling is strong at early times [12] . Such a strong- 
coupling scenario has previously been dismissed |12) . 
since the free-field behavior of electromagnetic waves 
would no longer be valid. We therefore consider here 
the alternative Lagrangian, 



d^x 



A'' J a + Co 



(10) 

in which the conformal factor l'^(4>) is moved outside 
the entire electromagnetic-sector Lagrangian. With this 
modification, the strong-coupling problem is alleviated. 
This Lagrangian could arise if is a dilaton field, al- 
though in that case we would not expect the conformal 
factor P{4)) to also multiply the mass term of the charged 
particle. 
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The effect of the coupUng function on Maxwell's equa- 
tions is straightforward. In the absence of charges or 
currents, only Ampere's equation is modified, to 



V X B 



1 d 



a^WE] , 



(11) 



where we have assumed W is solely a function of (con- 
formal) time. Faraday's law remains unaltered: 



V X S = - 



a? drj 



(12) 



A homogeneous magnetic field therefore scales with con- 
formal time as \B\ cx rf' , while a homogeneous electric 
field scales as \E\ cx ?y2+2»\ magnetic- and electric- 
field energy densities scale as pb — WB'^/Stt oc 
and = WE'^/Sn cx r]'^+^"'. Thus, for n = 2 and 
n ~ —3 (special cases we will consider below), the 
magnetic-field energy density remains constant or decays, 
respectively. The electric-field energy density decays for 
n — 2, but it grows, as pE cx 77^^, for n = —3. In 
this latter case, the energy density in the electric-field 
component of the quantum-mechanically induced elec- 
tromagnetic fields will, if infiation goes on long enough, 
ultimately dominate the energy density, ~ 3i?|/(87rG), 
in the inflaton. As we will see below (see also Ref. |12|). 
this then severely restricts the number of e-foldings of 
inflation. We will thus ultimately discard the n = —3 
case. 



C. Quantum Fluctuations of the Magnetic Field 

The action for the free field theory is 
1 



= I d^Sx [mf Q^f - \{d,A, d,A,f^ , 

(13) 

in the Coulomb gauge, where Ai is the vector potential. 
The Latin indices here are contracted using the spatial 
part of the Minkowski metric. Defining the vector field 
Vi = lijfjAi we can bring the kinetic term to canonical 
form, whereby 



Sera = i drj d^x 



1 



V 



-V. 



(14) 



after some integrations by parts. The quantized field Vi 
is expanded in terms of time-dependent mode functions 



d^k 

(27r)3 



e*"^^fe(ry)ei"Hfc)6.(fc)+h.c. 



(15) 



where 6, 6^ are annihilation and creation operators sat- 
isfying [ha(k),hi,(k)] = {2nf5,,a,5D{k-k'), where e^^^ 
is the polarization vector, a sums over the two linear- 
polarization states, and ^„ ef\k)e^"\k) — 6ij — kikj 
which further ensures transversality as a consequence of 
the gauge choice. Canonical quantization means that the 
vector field and its conjugate momentum V.^' satisfy the 
commutation relation. 



[Vi{x,T]),V-{y,'q)] = id^jd{x-y), 



(16) 



which results in a constraint to the two linearly inde- 
pendent solutions to the mode equation. The scalar field 
contributes an effective time-dependent mass term to the 
vector field, so that the mode functions obey the equa- 
tion, 



k'- 



1'' 



Vk = 0, 



(17) 



where /"// = n{n + 1)/??^, is positive for n > or n' ~ 
n+1 < 0. At high frequencies, k\ri\ ^ 1, the solutions are 
oscillatory, but at low frequencies the scalar field causes 
solutions to grow as Vk cx jTyj^^". The normalized 

solution, having positive frequency in the remote past, 
rj — > —00, for k\r]\ ^ 1, is 



2 V2k 5+"^ " 



(18) 



where Hn{x) is a Hankel function. In this case, the two- 
point correlation function for the magnetic field is 



{B{x,fj) ■ B{y,rj)) = 



1 



0(77)"* \y dx^dy^ dx^dy^ 
X (A,(x,r/)Aj(j/,r/)) 

^e-(-^V,(fc), (19) 



where 



PB{.k) 



2"P" 



2n 



{-k^)'H^^lj-k^)HYlJ-krj), 

m 

is the magnetic-field power spectrum. In the unamplified 
case, corresponding to n = 0, we have P^B^ = k{HirjiY at 
the end of inflation; the correlations in this case are then 
the usual vacuum-fluctuation correlations. Production of 
classical long-wavelength magnetic flelds occurs for n > 
or for n' = n + 1 < 0. To treat both cases with a single 
expression, we define ub = 4 — 2n for the case n > and 
ub — 4: + 2n' for n' = n-|-l < 0. Consequently, the power 
spectrum is 



(2) 



Pb 



r(^)2 



-kmnr' 



P. 



(0) 



B ' 



(21) 



for modes outside the horizon at the end of infiation. 
Since fc|?7/| <^ 1 for modes outside the horizon at the end 
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of inflation, and since ub — 4 < 0, the amplified ratio 
Pb/Pb^ can grow quite large. 

The mean-squared magnetic-field power in long- 
wavelength modes at the end of inflation, per logarithmic 
interval, is 



where the flnal term in the top line ultimately vanishes 
due to the transversality condition k ■ e. The energy- 
density power spectrum consists of two terms, a kinetic 
term due to the electric field and a spatial-gradient term 
due to the magnetic field. 



din fc 



hb 



Hji-k^i/2r^. (22) 



A scale-free spectrum ub — can be achieved for 
72 = 2, -3. Using (Gauss)V87r = 1.91 x IQ-^o QgV^ 
Mpc = 1.56 X 1038GeV"\ estimating \r]i\ ~ lO'^^ Gpc 
(consistent with Hj ~ 10" GeV and z/ ~ lO^^ for the 
redshift to the end of inflation), and then redshifting to 
the present day, we find 



dlnfc 



10 



-18-24. 3rii5 



r( 



5 — n£ 
2 



r(5/2)2 VMpc 



G^. 

(23) 

li riB = or n = 2 or —3, then the field strength is 
roughly 10^^ G, which may be sufficient to explain the 
observed astrophysical and cosmological magnetic fields 



D. Energy Density of the Magnetic and Electric 
Fluctuations 



The same magnetic-field spectrum is obtained for two 
values of the index n. However, the time evolution of the 
coupling function 1(77) breaks the usual duality between 
electric and magnetic fields, and the electric-field energy 
density may in some cases increase, as discussed above. 
We require the energy density in superhorizon modes of 
the electromagnetic fields to be smaller than the energy 
density of the inflaton, and thereby derive now a restric- 
tion on the allowed values of n and Hj. 

The stress-energy tensor that appears as a source for 
the Einstein equations is 



I\<f>) (^g^pF^^^F^f - Igi^'^F^pF'^P^ . (24) 



The energy density observed in the cosmic rest frame is 

(Pk 



PEB = ^^^{A\A\ + {d,A,){d,A,)~{d.A,){d,A,)) 



x^y J (2TTy 



Ppik), 



(25) 



Ppik) 



3(277)3 

2a4(77) 



(26) 

The integral over wavenumbers runs from fcmin = ~^lvs 
to fcmax = ^1/^77 where 775 = rjje^' is the conformal time 
at the beginning of Nj e-foldings of inflation, thereby 
spanning the range of wavelengths that have exited the 
horizon by the time 77. The pattern of behavior distin- 
guishes two regimes. 



PEB ^ Hf X 




2{\n\-l) 



for 
for 



\n\<2, 
\n\ > 2. 



(27) 

In the flrst case, which includes the scale-free solution 
n — 2, the energy density is simply proportional to Hf 
which is always subdominant to the inflaton energy den- 
sity. However, the second case, which includes the other 
scale-free solution n = —3, places severe restrictions. 



H<2+— In^, 
' ' Ni Hi' 



(28) 



on the index 



Since observational constraints limit 



Hj < lO^^A/p, then to achieve at least 60 e-foldings 
of inflation, the index is bounded by |n| < 2.2, thereby 
eliminating the case n = —3. At the value n = —2.2, 
Eq. ( 23 1 tells us that the magnetic-field strength on Mpc 
scales is roughly 10~'^'^G. The case n = 2, however, safely 
satisfies the above bound and yields a nG magnetic field 
as we have shown. 



III. CORRELATION OF MAGNETIC FIELDS 
AND SCALAR FLUCTUATIONS 

We now evaluate the {5(f>)BB correlation making use of 
the in-in formalism |29| . After splitting the Hamiltonian 
into a free part plus an interaction part i?int, we may 
evaluate, to first order in perturbation theory. 



(^j^i^^v)Ai{y,v)Ajiz,Tj) 



drji 2 Im 



(50 

^^int im ) v)At {y, v)Aj [z, rj) 



The interaction Hamiltonian is 



-2n 



{A^'~-{d,A,^d,A,y 



(29) 



(30) 



Using Eqs. Q and (15), we find that the expectation value on the right hand side of Eq. (29) is 



Hint{tl') Vi)My, 'ni)Aj{z, f]i) 



The functions K^]'' and K^f^ are defined as 



(Pki d?k2 d^h 
(27r)3 (27r)3 (27r)3 



-2n 



(32) 
(33) 



where we indicate the scalar mode functions of the vector potential as Akii]) = Vk{r/)/ I{r]). Plugging Eqs. (32|-(33) 
into Eq. (31 1, we find 



^^{x,7jj)AMvMjiz,Vi)'^ - /n ^e^^~'^-^+~'^-y+^'^-'H2nmki + k^ + hW^j, (34) 



M 



(35) 



where we introduce the integrals 

/•OO 

Ii = Im / d/iui(i + ^ui)(-i + Ui)e"'i(^-iV"^" 



dpi 
Im 



r(2) 



r(2) 



(36) 



(37) 



While I2 and the magnetic-field power spectrum are both invariant under n I + n, Xi is not. This is not surprising 
since the interaction Hamiltonian is not invariant under this operation. In the above, we have defined fj, = 77/777 and 
Ui = —kiiji for i = 1, 2, 3. 

The three-point correlation function for the scalar field with the magnetic field is obtained from 



— {x,r)i)B{y,rii) ■ B{z,rii)) = -2Im / drj (Hintiv) j^ix,ili)B{y,rii) ■ B{z,rii) 



1 



92 



0(777)4 y '^dy'^dz'' dyidz' J \ M 



(38) 
(39) 



After some calculations, the final result is 

f^ix,rji)B{y,r^j)-B{z,m))- J ^^^e^^~^-^'+~^^-y+~^^^^^^ (40) 



where 



^3(^1,^2,^3) = 



^2 fHi 



I k2h 



8 \M J a{r]iY kf 



^(2k2-hli + {l + {k2-hf)l2 



(41) 



From statistical isotropy, the function -P3(fci, fc2, ^3) de- and we have used fc2 • ^3 — (ki ^ ^2 — fc|)/(2fc2A:3) in 
pends only on the magnitudes of the three wavevectors, Eq. (41). Eqs. (40) and (41 1, along with Eqs. (36)-(37l, 



6 



form the main results on which our subsequent analysis 
is based. 



IV. ANALYSIS OF CROSS-CORRELATION 

We would like to analyze the cross-correlation between 
the primordial magnetic field and the scalar field to de- 
termine if there is any imprint or unique signature that 
would indicate the scalar field's role in the amplification. 



In the case of most interest, n = 2, the integrals yield 

36 Ml / o , ,N 

^1 - 3 3 2 - U1U2U3 - U}{uiU2 + U1U3 + U2U3)) , 

IT U'2U^(jJ 



I7 



(48) 



(49) 



A. The Amplified Cross-Correlation Power 
Spectrum 

To start, we calculate the cross-correlation power spec- 
trum for several trial cases. The integrals Ii^2 can be 
evaluated analytically for integer values of n. In most 
cases, the results are cumbersome, so we assume Ui <^ 1 
after carrying out the integrals in order to shorten the 
expressions. For example, for n — 0, 



4 ui{ui + uj) 



(42) 



where w = ui + 1*2 + U3. Plugging in these results, we 
find 



^3(^1,^2,^3) 



{Hi/Mf {2ki +k2 + ksKh - fc2 - 



8kik2k3 



(43) 

where we have used cos — k2-ky, — (fc^ — A;|)/(2fc2fc3) 
for the angle between the vectors /c2 and ^3. 

Amplification occurs for n > and n < —1, so that for 
comparison we consider integer cases n = I and n — —2 
whereupon the integration simplifies. For n = 1, we find 



I9 



4 ui(ui + uj) 



(44) 



4 Ui{uJ^ + u\uJ — UiLO^ — U2U3UJ — U1U2U3) 



(45) 



For n = -2, 

+3m^(U2W3 -I- U2U3 — U!'^) — U2U^{UJ + Ui)) ,(46) 
4 Ui{ll>^ + ufiU — U\ijP' — U2U3UJ — U1U2U3) 



4t o"o „ (3uluj'^ (7 + In cj) + 3ulu2U3 



(47) 



Since |fcr//| <C 1 we have discarded subdominant terms 
from the above results. We note that the In a; term above 
results in a large numerical coefficient, since — fcr?/ 
lO"^'^ for modes that are just entering the horizon to- 
day. These expressions are inserted into Eq. (41) to find 
the cross-correlation power spectrum. 



B. The Real-Space Cross-Correlation Coefficient 



Our next step is to determine the dimensionless magni- 
tude of the cross-correlation; i.e., how strongly does the 
magnetic-field energy density correlate with the scalar- 
field perturbation? We thus now calculate the zero-lag 
cross-correlation ((J0)i?^) in units of ((50)^)^^^ (^^)- 

This cross-correlation amplitude, evaluated in the co- 
incidence limit, can be evaluated as follows. Starting 



eliminate the 



de 



from Eqs. (40|-(41), we evaluate the ki integration to 

The remaining integrand 



ta function. 

depends only on the magnitudes fc2, k^, and 0, the angle 
between the two vectors: 



A'f /■ 

{{54>)B^) = ^J kldk2kldk3d{cOSe)P3{k^,k2M 

(50) 

where fci = (fc| -I- fc| -I- 2^2^:3 cos 0)^/^. However, we can 
replace the 9 integral by fci, whereby 



{{5<t>)B 



M 



k2dk2 k^dk 



|fc2-/C3l 



fclrffcl -P3 (^1,^2,^3 )• 
(51) 



Since the integrand is invariant under the exchange of ^2 
and ^3, we can replace P3 — 2P39{k2 — k^) and remove 
the absolute-value sign from the lower limit of integra- 
tion. We implement cutoffs at both large and small fc, 
for the ultraviolet and infrared divergences that arise in 
both the scalar and magnetic-field spectra. The cross- 
correlation for 71 = and n — 2 are 
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((S6)B') ^ X (1^^ - fi) ' - = 0' (52) 

^' ' - \&Tr^a^{m)\M ) |7774(l00 + 241n3r-721n2rln(-fc„,,,77,)), n = 2, ^ ^ 

where r = fcmax/^min, and fcmax and fcmin are upper and lower bounds on the run of wavevectors. In practice, we 
expect to Unk the minimum wavevector with the Hubble scale, /cmin — 27riJo, and the maximum wavevector with 
some galactic scale, fcmax — 27r/A where A ^kpc. Since 1^77/ 1 1, we have discarded subdominant terms from the 
above results. The dimensionless cross-correlation coefficient X^^b'^i formed from the ratio of the cross-correlation 
with the root-mean-square amplitudes of the scalar and magnetic fields gives 

'^"^ "(^0US?m."U(|f)(25 + 61nV-181nVln(-fc,„a.^,))Mn3r, n = 2. ^"'^ 



Considering a sufficiently wide range of scales, e.g. r > 
10*, then X{n = 0) ~ {Hi /M)^/\Kr /-k and X{n = 2) ~ 

8(i7//A/)\/b7ln(-fcmax??/)/7r. Using -fcmax??/ ^ 10"^^ 
then the cross-correlation coefficient in the presence of 
the amplification mechanism is enhanced by a factor 
of ^ 500 over the case without the magnetic-field am- 
plification mechanism. When the full range of infla- 
tionary length scales is taken, r ~ 10^^, then X(n = 
2) ~ 2 X 10'^(i?//M). Since the cross-correlation coef- 
ficient cannot exceed unity, we infer an upper bound of 
Hj/M < 5 X 10""* which is consistent with naive expec- 
tations based on an inflationary scenario. 



scalar-field and magnetic-field power spectra are now 

((^0/M)2) = ^e™-(-~-^~)/^P5^, (57) 



= V-^Ps^/M^, (58) 
(B^) - 5]e™-(-~-^^/^,5a,,,,,P^, (59) 



V-^Pe 



(60) 



so that Pstj, is dimensionless and Pb has units of 
(energy)"*. The three-point function becomes 



C. The Behavior in Fourier-Space 

We now evaluate the triangle-shape dependence of the 
full three-point correlation function in Fourier space. To 
do so, we evaluate a ratio of the form. 



P3(fci,fc2,fc3) 

y/Psc^iki)PB{k2)PB{hy 



(54) 



to normalize the cross-correlation power spectrum. How- 
ever, since this ratio is not dimensionless, given our 
Fourier conventions, we go to a discretized Fourier trans- 
form, 



d^k 1 



(27r)3 V ^ 



(55) 



and likewise replacing the Dirac delta function with a 
Kronecker delta. 



M 



B') = 



^ ^i{ni-x-\-n2-y+n^-z)/Lp^ 



ni+n2+^3— 



Pa - V-^P:i, 



(61) 



where P3 has units of (energy)"*. We can now build a 
dimensionless cross-correlation coefficient. 



Cn = 



P3ini,n2,n3) 



Ps,p[ni)PB{n2)PB{n3) 



(62) 



where rii for i — 1, 2,3 are the magnitudes of vectors ni 
that form a closed triangle. 

For isosceles triangles with n2 — n^, the correlation C„ 
obtained for the case n = and ?i = 1 is 



Co = 



1 Hi {ni + n2){ni - 2n2) 



87r3/2 M 
1 Hi 



3/2 3 
N 



167r3/2 M (n, + 2n2)n?/'ni' 
N ^ nl + 2n\n2 - 2n\nl - Qn\nl 



(63) 



{2nf5{ki + k2)^V5n,,n^. 



(56) 



We presume a maximum length, L, so that the volume 
is y = and mode numbers are ki = 2'i:ni/L. The 



where — 1 < cos 6 = i^-T- — 1 < 1. An expression for C2 
^ "2 

is easily calculated, but the result is rather long and un- 
enlightening. The behavior of C„(cos6') for n = 0, 1, 2, 
and —2 is illustrated in Fig. [l] 
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FIG. 1: The ratio C,i(cos S)/C„(— 1) is shown for n = 0, 1, 2, and —2, as functions of cos^. In the n — —2 panel, the dashed 
line indicates where the absolute value has been taken. In the n = 2, —2 cases we have used 2Tvni\rii / L\ ^ 10~^^ corresponding 
approximately to a Gpc length scale. Note that the case of cosmological interest is n = 2. 



We find that there are two interesting limits for isosce- 
les triangles with n2 ~ n^, first a squeezed triangle, with 
1 < ni <C ^2 or ~ vr, and second a flattened triangle, 
with 712 = ni/2 OT 6 — 0. For the squeezed triangle we 
find the universal result. 



C„(cos7r) = 



V2 Hi 
(27rni)3/2 M ' 



(64) 



for all values of n, as borne out by numerical integration 
for non-integer n. We suspect that this triangle con- 
figuration, with small ni and large n2, 71.3, dominates the 
integration in Eq. ( 53 ) , as a way to help explain the simi- 



larities seen in the real-space cross-correlation coefficients 
for different values of the index n. 



The result, Eq. (64), suggests a natural reference point. 



so that a general expression for the discretized Fourier- 



space dependence of the cross-correlation is 



C„ (cos 6')/C„ (cos tt) — 



TT 2cos6iIi + (l 



cos^ e)i2 



4eni I 



(1) 
l/2^ 



,(^"3)1 

(65) 

where e = —2TTrij/L <^ 1. 

For a flattened triangle, we have Cq (cos 0) / C„ (cos tt) — 
0, Ci(cos0)/C„(cos7r) = 3, and C_2(cos0)/C„(cos7r) — 
12 (2 — 7 — ln(2e)), where 7 is the Euler-Mascheroni con- 
stant. Note that the cross-correlation vanishes for the un- 
amplified case (n = 0), but grows large for n = —2, where 
the argument of the log is ^ 10^^^ for modes entering the 
horizon today. The behavior of C„(cosO)/C„(cos7r) as a 
function of n is shown in Fig. [2j 

To show the full Fourier-space triangle dependence of 
the cross-correlation, we define the quantity 



R = 



n2 \ Cn (cos 9) 
ns J C„(cos7r) 



(66) 
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FIG. 2: The ratio C„(l)/C(— 1), the ratio of the discretized Fourier-space cross-correlation coefBcients for the flattened triangle 
to that of the universal result for the squeezed triangle, is shown as a function of n. No amplification, n = 0, yields zero 
cross-correlation. Hence, the flattened triangle may be used as an indicator of an ampliflcation mechanism. The ratio is 
negative along the dashed line, where we have taken the absolute value. We have set 2iTni\rii / L\ ~ 10~® for ease of numerical 
computation; using 27rni|r;j/L| ~ 10~^^ to represent Gpc scales boosts the curve up to 10'* near n = ±2. Note that the case 
of cosmological interest is n = 2. 



and introduce the variables 2:23 = n2/n^ and X13 = 
ni/n^, where < X23 < 1 and 1 — 0:23 < 2:13 < 1 + 2:23 
covers the full set of triangles. The behavior for the cases 
n = 0, 1, 2, and —2 are shown in Fig.[3| The Figure helps 
illustrate the difference between the amplified and unam- 
plified (n = 0) cases, and also shows that the maximum 
value of R in the amplified case occurs for the flattened 
triangles, corresponding to the line ccis = 1-1- 2:23, along 
which 9 — Q. Squeezed triangles, where — n, are lo- 
cated along 2:13 = 1 — 2:23. 



The amplitude of the zero-lag, real-space cross- 
correlation is dominated by the cut-offs to the range of 
Fourier modes. At finite Fourier mode, however, the 
flattened isosceles triangle produces the largest cross- 
correlation amplitude and hence the best opportunity to 
observe the signature of the mechanism responsible for 
amplifying the primordial magnetic field. 



V. CONCLUSIONS 

We have considered a toy model in which a scalar field 
is coupled to electromagnetism in a fixed de Sitter back- 
ground. The homogeneous time evolution of the scalar 
field breaks the conformal invariance of electromagnetism 
resulting in quantum production of magnetic fields in 
addition to quantum production of scalar-field fiuctua- 
tions. We then calculated the cross-correlation between 
the scalar field and the magnetic-field energy density. 
The dimensionless cross-correlation coefficient is propor- 
tional to the ratio Hj/M, which must be small if the 
effect of the scalar-field perturbation on the electromag- 
netic part of the Lagrangian can be considered small. 
However, this small quantity may be multiplied by a nu- 
merically large (< 500) coefficient suggesting a possibly 
strong (even order-unity) cross-correlation. 

We also studied the full triangle-shape dependence of 
the three-point correlation function in Fourier space. We 
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FIG. 3: The quantity R, defined in the text as the ratio of the Fourier-space cross-correlation coefficient to that of the 
universal result for the squeezed triangle, times a factor 3:23, is shown as a function of the triangle side lengths. We have set 
2Tvn3\rii / L\ ~ 10^^^ for the cases n = ±2. Note that the case of cosmological interest is n = 2. 



find that it is nonzero for squeezed triangles (wherein 
the short Fourier component is that associated with the 
scalar-field mode), but may be considerably larger for 
flattened triangles (where the long Fourier mode is as- 
sociated with the scalar field and twice as long as those 
associated with the magnetic field). These shape depen- 
dences may be useful if such correlations are to be sought 
in the data. 

Although we treat it as a toy model, our calculation 
provides the correlation between the curvaton and mag- 
netic fields if the scalar field is identified as the curvaton. 
If primordial perturbations are further due to curvaton 
fiuctuations, then the scalar-field-magnetic-field cross- 
correlations derived here describe density-perturbation- 
magnetic-field correlations in the Universe today. If the 
scalar field is the inflaton, then there are additional steps 
to relate the scalar-field perturbation to the density- 
perturbation amplitude in the Universe today [36]. 

The cross-correlation between primordial-seeded den- 
sity perturbations and magnetic fields may be amenable 
to detection through the cosmic microwave background 
(CMB). Cosmic magnetic fields present during the re- 
combination era contribute to the CMB temperature and 
polarization signals. (See Ref. [TS] for a detailed study.) 



Magnetic fields along the line of sight further distort the 
CMB by converting E-mode polarization into B-niode po- 
larization, through Faraday rotation [TBti20j . Primordial 
magnetic fields may also leave a non-Gaussian imprint on 
the statistics of the anisotropy pattern [21- 25J. Current 
observations set the upper bound on a primordial mag- 
netic field at the nG level (We also note that 
there have been claims of a lower bound on an extra- 
galactic field [SO].) 

The correlation may also be accessible through a com- 
bined survey of large-scale structure and Faraday rota- 
tion [31]. The proposed SKA telescope, which is pro- 
jected to be sensitive to variations of 0.1 nG across 100 
Mpc, and LOFAR, which aims to explore the nG fields in 
intergalactic media |32H35j , may offer more direct means 
to probe for the cross correlation. The detectability of 
the effect studied here will be the subject of future work 
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